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For systems exhibiting a second-order phase transition 
with a spontaneously broken continuous 0(A r )-symmetry at 
low temperature, we give a criterion for judging at which tem- 
perature Tk long-range directional fluctuations of the order 
field destroy the order when approaching the critical temper- 
ature from below. The temperature Tk lies always signifi- 
cantly below the famous Ginzburg temperature Tc at which 
size fluctuations of finite range in the order field become im- 
portant. 

1. Although the fluctuation behavior of systems in a 
second-order phase transition is universal in the imme- 
diate vicinity of the transition, and the theory is well 
established jij, there are still disputes concerning the 
dominant fluctuation mechanism which drives a num- 
ber of important phase transitions. A typical example is 
the question whether the superfluid transition in liquid 
helium is initiated by the proliferation of vortex lines, 
which are the defects in all systems with pure angular 
fluctuations of a complex order parameter, or by the size 
fluctuations of the associated complex order field (f>(x). 
At first sight, this question may seem meaningless since 
all critical exponents of the transition can be calculated 
with great accuracy either from an order field theory 
with (^-interaction ||, where size and directional fluc- 
tuations seem to be equally important, or from Heisen- 
berg model on a lattice, which contains only directional 
fluctuations. The equivalence can easily be proved for su- 
perfluids, where the Heisenberg model reduces to an XY- 
model. After a duality transformation, the XY-model 
can be reexpressed as a sum over a grand-canonical en- 
semble of non-selfbacktracking vortex lines ||[|] whose 
proliferation completely describes all properties of the 
superfluid transition. They produce the same critical ex- 
ponents as a complex </> 4 -theory, and the reason for this is 
simple: the XY-model may be converted into a complex 
(/> 4 -theory by a transformation of integration variables in 
the functional integral of the partition function Q . 

The question whether angular or size fluctuations drive 
a phase transition does therefore not concern the imme- 
diate vicinity of the transition where the properties are 
universally governed by a critical exponent w. ft is only 
a meaningful question in the precritical regime of the 
transition, where a mean-field description of a system 
breaks down. There it possesses an answer similar to 
the Ginzburg criterion which estimates the temperature 
range where this happens due to size fluctuations of the 
order field. The new criterion to be presented in this note 
will tell us where the breakdown is caused by directional 



fluctuations. In the case of a complex order field, these 
lead to an early proliferation of vortex lines before size 
fluctuations become large. 

The new criterion serves to understand the dominance 
of directional fluctuations in the recently-discovered 
restoration of continuous symmetry in Gross-Neveu || 
and Nambu-Jona-Lasinio models and, more impor- 
tantly, the generation of a pseudogap phase above the 
superfluid phase in strong-coupling superconductors . 

2. In order to lead up to the new criterion, we briefly re- 
call the relevant features of the field-theoretic approach 
to the critical exponents in the immediate vicinity of 
second-order phase transitions |Q. Since critical prop- 
erties are caused by the long-wavelength fluctuations of 
a system, it is sufficient to identify these, assign to each 
of them a real order field 4>a{x) : A = 1, . . . , N, the local 
generalization of Landau's order parameter [8j, and set 
up a Ginzburg-Landau energy density. Its fluctuations 
are studied with the help of a functional integral over all 
field configurations, weighted by a Boltzmann factor of 
the total field energy. 

At a moderate temperature distance \T — T c \ from the 
critical temperature T c , the functional integral may be 
evaluated by the saddle point approximation. This is the 
mean-field regime of the field theory. When approaching 
the critical point, the fluctuations of the order parame- 
ter increase. They become important in a temperature 
regime /STq around T c , whose width was first estimated 
by the famous Ginzburg criterion |9j: It is the regime 
where the fluctuations of the order parameter in pockets 
of coherence size reach into the normal phase. 

The Ginzburg-Landau field energy density is found 
from phenomenological studies of a system in some neigh- 
borhood of the critical temperature ||. One expands 
the energy density in powers of the order field and its 
derivatives. Then one identifies the temperature interval 
where all expansion terms are irrelevant except for those 
appearing in a simple </> 4 -theory. At some mean-field crit- 
ical temperature T C MF ^ T c , correlation lengths become 
infinite in the mean-field approximation. The mass term 



the bare square 
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goes through zero linearly in T / , i.e 
mass m 2 is proportional to r = T/T^ F 

As the system enters the critical regime, fluctuations of 
the order parameter become important and must be ac- 
counted for. This is done by perturbativc methods. Each 
physical quantity of interest is expanded into a Taylor se- 
ries of the reduced interaction strength g/m 4 ~ D , and a 
resummation of these in the limit m —> allows us to 
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extract power laws m powGr , which determine the critical 
exponents [^0| . 

We shall now demonstrate that in systems with a con- 
tinuous symmetry, the Ginzburg criterion is not sufficient 
to characterize the temperature range where fluctuations 
are important. A system will exhibit strong directional 
fluctuations before the size fluctuations of the order pa- 
rameter become noticeable, leading to a phase transition 
at a temperature Tk far below the Ginzburg temperature 
T G - 

3. For simplicity, we restrict the argument to N order 
fields <f>A with 0(iV)-symmetry in D dimensions. With 
a convenient choice of field and mass normalization, the 
Ginzburg-Landau energy density in D dimensions may 
be written as 

£ (^,a^)=^ IJ {aV[a^( a; )] 2 +T0i(x)+| [<t> 2 A {x)] 2 ) . 

From here on we use natural units with ksT^ = 1. 
The fields have zero engineering dimension, a denotes 
some microscopic length scale of the system, usually 
the size of atoms or molecules, and g is some inter- 
action strength. The parameter a specifies the zero- 
temperature coherence length of the system in units of a 
as being £o = aa/y2. This can vary greatly from system 
to system. In superconductors, for example, a it can lie 
anywhere between a few thousand, and less than ten in 
high-temperature superconductors. 

In this note, we shall only be concerned with the de- 
struction of the ordered state which lies below the crit- 
ical temperature where r < 0: There the fields in the 
energy density fluctuate around an ordered ground state 
with a constant vector (<fr A ) = &a = (<f>) Na = &Na 
in field space, whose direction vector Na breaks spon- 
taneously the 0(A r )-symmetry, and whose magnitude is 

^ = V ^\ = \/ ~ T 1 1 9i where the energy density is mini- 
mal, fluctuating around the condensation energy density 
£o = e(^A)O) = —r 2 /4:ga D . The temperature-dependent 
coherence length £ = aa/y2[r] describes the range of 
the size fluctuations of the order field. 

The magnitude is estimated by assuming the field to 
live in patches on a simple cubic lattice of spacing £j = l£, 
choosing eventually a spacing parameter I = 2 to ensure 
the independence of the patches. Then 

(foOr)-*] 2 )/* 2 = l 2 - D {2\r\) D ^ga- D v 1 {Q), (1) 

where v^ 2 (0) = j! w d D k/[£? =1 (2 - 2cosfc 4 ) + to 2 ] is 
the lattice Coulomb potential of reduced mass to. For 
D = 3, 4, uf(0) has the values [p| 

uf(0) w 0.1943, 0.1270, 1/2D. (2) 

Mean-field behavior breaks down if ([!]) is of the order 
unity, which happens at the reduced Ginzburg tempera- 
ture 

\t g \^[Kv?{0)/1] 2 '^- d \ D<4. (3) 



where 

K = 7Pt*-ig/oP, (4) 

i.e., at a Ginzburg temperature Tq = T C MF (1 — \tc\). 
Ginzburg, in his original paper H, estimated v®(0) in 
three dimensions by an integral J d 3 p/(2Tr) 3 (p 2 + 1) 
(27T 2 ) -1 f£ dpp 2 /(p 2 + 1) s» 1/47T, and assumed / = 1, 
which lead to \tq\ ~ (g/a 3 ) 2 /8ir 2 . In old-fashioned 
type-II superconductors, \tq\ can be as small as 1/10000 
p2| , which explains why conventional superconductors 
are well described by mean-field theory. In modern high- 
T c superconductors, on the other hand, Ginzburg's esti- 
mate leads to \tq\ w 0.01 [^3[ , such that critical expo- 
nents become observable. 

For D > 4, the right-hand side in (jjj) decreases when 
approaching the critical point, so only mean-field behav- 
ior is observed. If D = 4 — e lies only slightly below four, 
the right-hand side of (||) behaves like |r| e / 2 , implying a 
good mean-field description until |t| is extremely small. 

4. The derivation of the new criterion is based on the 
observation that the kinetic term defines a second, com- 
pletely independent, energy scale of the system. For its 
identification, we split the fields according to size and 
direction in O(N) field space as 4>a — 4> n A, n 2 A = 1. 
The directions ua describe the long-range fluctuations 
of the Goldstone modes. Sufficiently far from the criti- 
cal regime, we may neglect the gradient term of the size 
4>{x), and approximate the energy density by 

e(<p,dn A ) = ^ 5 {aW(aO [dn A {x)] 2 + T<t> 2 (x) + g -^{x)}. 

The fluctuations of the Goldstone modes are controlled 
by the gradient term whose magnitude depends on the 
size $ of <p at the minimum of the potential. The gradient 
energy density is 

e nA (dn A ) = ^_ 2 [dn A {x)} 2 , (5) 

with = l3{<P)=a 2 {Zi/a) D - 2 <S> 2 =a D l D - 2 /(2\T\) D / 2 - 2 g. 
This is the second energy scale. It measures how much 
energy is spent when reversing the direction vector ua 
over the distance £/, and is called the stiffness of the 
directional field. 

From studies of 0(iV)-symmetric classical Heisenberg 
models it is known that directional fluctuations disorder 
a system if the bending stiffness drops below a certain 
critical value j3 CI . For large N, this value can easily be 
estimated by a simple manipulation of the functional in- 
tegral for the partition function associated with the en- 
ergy (||). It may be written as frl4| 

Z n , = Jvn A V\e~ {P/2iD ~ 2) I dD *{[ dn A(x)] 2 +\(x)[n\(x)-i]} _ 

where the unit length of n A (x) is enforced by a Lagrange 
multiplier field X(x). Integrating out the n^aO-fields 
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leaves us with a pure A(a;)-field theory, with an energy 
functional 



E[X] = yTrlog[-a 2 + A(x)] 
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d D x\{x). (6) 



For large N, the fluctuations of the A-field are frozen, and 
the disordered state has an energy density given by the 
extremum of (^), where A(ai) is a constant satisfying the 
gap equation (3 = (3\ = Nv®(0). The order is destroyed 
at a critical stiffness at /3 cr = (3q. On a simple cubic 
lattice, we find in three, four, and large-_D dimensions 

& 

(3 cr = Nvg(0) ps N 0.2527, Y0.1549, N/2D, (7) 

respectively. Formula (Q) is reliable only for large N. 
However, Monte Carlo simulations show that the crit- 
ical value ^ can be trusted already for D = 3 and 
N = 2 to within about 10%, where simulations yield 
/3^ c ps 0.45 (see |l5|), in good agreement with the value 
0.5054 from (0). The simulations are done by putting 
the Heisenberg model on a lattice of unit spacing, so that 
the energy density for N = 2 takes the XY-model form 
£n A {dn A ) ps 3J2p,=i ...,d\ 1 ~ cos ^nl( x )]y where V 4 de- 
notes the lattice gradient in the ith coordinate direction, 
and 7 = arctanri2/ni. Since the quality of the approxi- 
mation increases with N and D, we can trust Eq. (Q) to 
within 10% for all N and D > 3. This accuracy will be 
sufficient for the criterion to be derived here. 
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FIG. 1. Softening of stiffness of XY-model derived from a 
self-consistent approximation a la Hartree-Fock-Bogoliubov. 
The dashed curve shows the mean-field approximation given 
by p 3 = a 2 /AD 2 (3 2 , fl = aI (a)/2Dh{a) [I n (a) = modified 
Bessel functions] which goes to zero linearly in |r|. The exact 
stiffness goes to zero like \T C — T^ 15-2 )", with the critical 
exponent v (ps 0.6705 for D = 3 Q). 

The critical stiffness can, incidentally, be also estimated 
by calculating its renormalized version from a sum of 
an infinite number of terms in a perturbation expan- 
sion. Expanding the cosine into a Taylor series, and 
calculating the harmonic expectation values of quartic, 
sextic, etc. terms, we find in a self-consistent approx- 
imation of the Hartree-Fock-Bogoliubov type that the 
stiffness has a renormalized value |l6| (3r = (3 c^ 1 ^ 2Di3r . 
This softens with increasing temperature 1//3, until (3 
reaches a critical value (3 CI — e/2D where /3r drops to 
zero (see Fig. |l|). In D = 3 dimensions, this happens at 
/3 cr = 0.4530 . . . , a value which is in excellent agreement 



with the Monte Carlo number /3„ C ~ 0.45. The predic- 
tion of such sharp drop is true only in two dimensions, 
as shown by Kosterlitz and Thouless [0. For D > 2 
it is an artefact of the approximations, and the exact p s 
goes to zero like \T C — T\( D ~ 2 ' V , with a critical exponent 
v ps 1/2 + (4 - D)/1Q + . . . . 

The estimate for the critical stiffness (Q) leads now di- 
rectly to the announced criterion: The phase fluctuations 
will disorder the system if the stiffness (3 in Eq. drops 
below the critical value (0), which happens at a reduced 
temperature 

\t k \ ps [NK v° (0) /lf/^- D) , D < 4, N > 2. (8) 

Thus we obtain the important result that 

\t k \ ps [iV^o(0)K(0)] 2 /( 4 - D )|r G |, £><4, N>2. (9) 

This implies that for all systems with N > 2, direc- 
tional fluctuations destroy the order before size fluctua- 
tions become large. They cause a phase transition below 
the Ginzburg temperature at Tk = T' C MF (1 — \tk\)- For 
D = 3, the relation becomes \t k \ ps 2.184 Y 2 |t g |. Thus, 
if the critical regime is approached in a 4 -theory with 
a well-formed mean-field regime, the transition is always 
initiated by directional fluctuations. In particular, the 
estimates for the critical regime of the high-|T c | super- 
conductors |llj] will receive a factor ps 9. 

The dominance of directional fluctuations is, of course, 
most prominent the limit of large N, and it is therefore 
not surprising that the critical exponents of the </> 4 -theory 
and the Heisenberg model have the same 1/iV-expansions 
in any dimension D > 2, as a pleasant demonstration of 
the universality of critical phenomena. 

By adding to the field energy density e(0, Bua) the 
energy density of directional fluctuations with the field- 
dependent stiffness (3 = (3(6) = a D 6 2 l D - 2 /{2\t\) d / 2 - 1 
we can study, as in Ref. ||, the combined energy den- 
sity in the disordered phase where the symmetry is re- 
stored but the average $ of the size of the order field 6 
is nonzero. 

5. How do we determine experimentally the fluctua- 
tion parameter K to estimate \tq\ and |tj<-|? In mag- 
netic systems, one measures the susceptibility tensor 
XAs{k) = J d D xe lkx ((j) J \(x)(j) B (0)} at wave vector k, and 
decomposes it into parallel and perpendicular parts as 
XA B (k) = ($A*B/$ 2 )x||(fc) + (Sab ~ <S> A <$>B/<S> 2 )xAk). 
The mean-field behavior of these quantities are \\\{k) « 
a 3 /(a 2 a 2 k 2 + 2\t\) and XJ-(^) ~ a 3 /a 2 a 2 k 2 . Combin- 
ing these at k — with the mean- field b ehavior of the 
spontaneous magnetization $ = y\r\/g, and with the 
temperature-dependent coherence length £, we see that 
the size of K can immediately be estimated from a plot, 
versus t = T/T c — 1, of either of the dimensionless exper- 
imental quantities 
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these being written down in in physical units. Note that t 
measures the temperature distance from the experimen- 
tal T c , in contrast to r = T/T C MF - 1. In the mean-field 
regime, where t ps r, K exp is constant, and can be in- 
serted into Eq. (||) to find the temperature Tk where 
directional fluctuations destroy the order. 

In superfluid helium we may plot, in analogy to the 
transverse susceptibility expression for K eKpi the quan- 
tity K cxp ps \t\ 2 - D / 2 M 2 k B T/£ D - 2 h 2 p s , where M is the 
atomic mass and p s the superfluid mass density, which 
at the mean-field level is defined by writing the gradi- 
ent energy (|) as (p s /2k B T)(h 2 / M 2 )[dn A {x)] 2 . In the 
critical regime, the three expressions for K cxp go uni- 
versally to zero like \t\ 2 ~ D ^ 2 , since £ cx \t\~ u , X||(0) pa 
\t\<*~ 2 >, k 2 X± (k)\ k ^ « \t\*», $ 2 ps \t\»( D - 2 +*\p s « 
\t\ ( -°- 2 >, with r\ ps [(TV + 2)/2(JV + 8) 2 ](4 - £) 2 + . . . . 

Experimentally, the superfluid density of helium for 
D = 3 shows no mean-field behavior a la Ginzburg- 
Landau down to T ps T c /4, such that the above for- 
mulas cannot properly be applied. Let us neverthe- 
less estimate orders of magnitude of a would-be mean- 
field behavior: p s /p ps 2|t| p^] , where p = Af/a 3 is 
the total mass density, with a ps 3.59 A Then 
the factor ksT c at T c = 2.18 K can be expressed as 
k B T c ps 2. 35^ 2 /Ma 3 With £ ps 2 A, we obtain 

an estimate K ps 1.2 a/£o ps 2. Inserting this into Eq. (||) 
and relation (^|), we obtain 



1/4, 



0.03. 
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The large size of \tk\ reflects the bad quality of a mean- 
field description. 
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